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Abstract 



Four new variants of the Computational Order of Convergence (COC) of a one-point it- 
erative method with memory for solving nonlinear equations are presented. Furthermore, the 

■ way to approximate the new variants to the local order of convergence is analyzed. Three of 
, the new definitions given here do not involve the unknown root. Numerical experiments using 

adaptive arithmetic with multiple precision and a stopping criteria are implemented without 
using any known root. 
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^ ! 1 Introduction 

'NT 

■ One-point iterative methods with memory for solving a nonhnear equation f{x) = 0, where 
O . / : / C M ^ M, and / is a neighborhood of the root o?, usually consider a sequence {3^n}ngN) 
^ I defined by 

Xn-i-l = 4>{Xn;Xn-l,---,Xn-j) , n > 0, (1) 



where (p is the iteration function. A sequence {xn} is said to converge to a with local order of 
convergence p G R, p > 1, if there exists the following limit 



"^I T log|e„+i| 

p = lim 



n^oo log|e„| 

where = — a is the error in the A;th iterate (see [Ill2])- This limit p is also equal to i?-order 
defined in For one-point method with memory ([TJ the error equation is: 

e„+i = Ce„^(l + 0(e-)), (3) 

where C is a real number, < cr < 1, and we will consider p > {1 + ^/5)/2. The nonzero constant 
C is called the asymptotic error. The local order of convergence of an iterative method in a 
neighborhood of a root is the order of the corresponding sequence. If it is p, then the method 
approximately multiplies by p the number of correct decimals after each iteration. That is, from 
jl]) we get logio kn+i| ~ P logio |e„|, for n large enough. 

In the next sections the way to approximate four new variants of the local order of convergence 
is analyzed and numerical experiments using adaptive arithmetic with multiple precision and a 
stopping criteria are implemented without using any known root for three of the four techniques. 
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2 Definitions and first result 



Next, we give the definitions of Computational Local Order of Convergence (CLOC) that is a 
variant of COC ([4j, 2000), Approximated Local Computational Order of Convergence (ACLOC), 
Extrapolated Local Computational Order of Convergence (ECLOC) and Petkovic Local Compu- 
tational Order of Convergence (PCLOC). These three last concepts are variants of ACOC, ECOC 
[7j and PCOC [9] respectively. After the work of Weerakoon and Fernando |4], many other au- 
thors have considered the COC in their research (see tSj-l2j and references therein). In all those 
papers the COC is used to test numerically the order of convergence of the methods presented. 
Considering ([2|) we provide a new parameter with lower cost than COC: 

Definition 1 The computational local order of convergence (CLOC) of a sequence {xn}n>o is 
defined by 

A„ = ^i^, (4) 

log I en- 1 1 

where Xn-i and x„ are two consecutive iterations near the root a and Cn = Xn — ol. 

Notice that the last definition has lower cost because we use the logarithm function applied to 
only one variable, say e„, instead of a quotient such as e„/e„_i which is used in [4J. 

The main drawback of COC and CLOC is that they involve the exact root a, which in a real 
situation it is not known a priori. To avoid this, we introduce three variants of CLOC that do not 
use the exact root. Firstly, we give a new parameter considering three consecutive points: 

Definition 2. The approximated computational local order of convergence (ACLOC) of a sequence 
{xn}n>o is defined by 

A. = ^^^, (5) 

log I Cn-ll 

where Cn = Xn — a^n-i • 

Secondly, in order to avoid the requirement of the knowledge of the exact root a, we consider 
three consecutive iterates Xn, Xn-i, Xn-2, and using Aitken's extrapolation we give the following 
approximation of a 2 

an = Xn """^ , n > 2, (6) 

where S is the forward difference operator, dx^ = x^+i — x^ and ([6]) is the J^-Aitken procedure [8]. 
Then, we can define a new approximation for the error e„ = x„ — 5^ and a new computational 
order of convergence: 

Definition 3. The extrapolated computational local order of convergence (ECLOC) of a sequence 
{xn}n>o is defined by 

A. = -i^, (7) 
log I e„_i| 

where Cn = Xn — ctn o^nd an is given by 

Finally, another way to avoid formulae involving the exact root a consists in using the values 
of two consecutive iterates. That is, from f{xn) and f{xn-i) the new computational order of 
convergence is: 

Definition 4 The Petkovic computational local order of convergence (PCLOC) of a sequence 
{xn}n>o is defined by 

X _ log I f{Xn)\ 

log I /(x„_i)| 
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This last parameter PCLOC is defined in honor of Petkovic who in [9l [10] consider, in analogy of 

Pn 



COC, the following value log |/(x.+i)//(x„)| 



log|/(3;„)//(x„_i)| ■ 

As we show below, for all sequence converging to a, with starting points x^j, . . . ,x_i,xo 

close enough to a, the values of A„, A„, A„ and A„ converge to p, when n — )• oo. 

There exist numerical problems where a huge number of significant digits of the solution is needed. 
Such problems require the use of methods with a high order of convergence together with adequate 
arithmetics. We compute with a multiple precision arithmetic or symbolic manipulators, as Maple, 
that allow us to work with an adaptive arithmetic, that is, to update the length of the mantissa 
in each iteration by means of the formula 

Digits := [p X (-log|e„| + 2)] , (9) 

where p is the local order of convergence of the method and [x] denotes the integer part of x. 
Notice that the length of the mantissa is increased approximately by the order of convergence 
p. In our experience, in order to guarantee all the significant digits required we add 2 units to 
-log|en|. 

Our first aim is to express as a function of tn+i- In a first approximation, from ([3]), if we 
consider e„+i = C Cn then we get = C^^^'' el/^i- Substituting this result in the second term 
of the right side of ([3]) we obtain e„+i = C Cn {l + 0{ e^^^)) , and 

eP = C-'en+i{l + 0{e^i'[)). 
Therefore, expressing in terms of e^+i we can state the following proposition: 

Proposition 2.1 Considering true the hypothesis in ^ we have 

en = C-^/^e„^^,(l + 0(e:/^,)). (10) 

3 Computational Local Order of Convergence (CLOC) 

A relationship between An and p is derived. In fact, we prove that A„ converges to p when 
en-i — 0. That is A„ ~ p, in the sense that lim — = 1. 

n— >oo p 

Proposition 3.1 // A„ is the CLOC defined in ^ and p is the order of convergence, then 

' log|g| 
\ log |en-i 

where C is given in (0/. 
Proof. To prove (jlip we express A„ in terms of e^-i. Taking into account ([3]) we take 

logle^l = log|Ce^_i(l+ 0(e„^_i))| 

= plog|e„_i|+log|C|+0(e„"_i). (12) 

Then 

log |e„| plog \en-i\ + log \C\ + 0{e^_^) 



A„ = p[ 1+ 0(_i^sM ) 1 (1^) 



log|e„_i| log|e„_i| 
and we obtain the assertion of the proposition. □ 

Notice that for the calculus of the CLOC ([!]) and for updating the adaptive arithmetic process ^ 
it is necessary to know the exact root a. In this case the following stopping criterion is applied: 

|e„| = |x„-a| < 10^", (13) 

where rj is the maximum number of correct decimals and 10"'' is the required accuracy. 
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4 Approximated Computational Local Order of Convergence (ACLOC) 

A relationship between and p is obtained. A new technique to update the number of significant 
digits in an adaptive multi-precision arithmetic is given and a new stopping criterion is suggested. 



Proposition 4.1 Let A„ he the ACLOC defined in ([^ and p the order of convergence, then 

X„., (1+0(^5^)), (14) 

V \pl0g\en-2\J J 

where C is given in 

Proof. From the following expression: 

log|en| = log|en-en-i| = log|en-i| +0(|e„/e„_i|), (15) 

and ([5]), we get 

- _ log|e„| log I e„_i| + 0(|e„/e„_i|) 



log I e„_i| log I e„_2| + 0(|e„_i/e„_2|) ' 
Applying ()12p we obtain 

^ ^ plog|en-2| +log|C| +0(e^_2) 
log I e„_2| + 0(|e„_i/e„_2|) 

where v = mm{p{p — 1),<7} and the proof is complete. □ 

Observe that from pT]) and the expressions of An and A„+i are identical. That is, if we 
approximate the theoretical value of the local order p by the computational values An and A^: 

P = An lb AAn, p = An+1 ± AAn+1, 

then AAn+i ~ AA„. Our numerical experiments confirm this relation. 
Repeating ([3]) twice we obtain 

e„ = C''+ie^'_2(l+ 0(e-_2)), 
and now, we write Cn/en-i in terms of en_2: 



n-2 ) 



Cn-l Cn-l - en-2 C e^_2 + 0(e;^^2 ) " en-2 

= Cer2(l + 0(en^-2)), (16) 
where r = min {p — 1, cr}. Moreover, we get 

en-2 = ^-(^+1)/"' e'Jp" (1 + ©(en-/")), (17) 

since, from m, en-2 = C-^P e^L^^ (l + 0{e^['[)) and en-i = C'^/p eU'' (l + 0(e"J^')). 
Substituting (fT7|) in (fTUj) we have the following proposition: 

Proposition 4.2 // we set Cn = Xn — ct and Cn = Xn — Xn-i, then 

en-CV(W)(^) , (18) 

where p is the order of convergence and C is given in 
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The result (jlSp allows us to substitute the error term in ([9|) by an expression that does not involve 
the exact root. Indeed, we implement the following adaptive multi-precision arithmetic scheme: 

n3 



Digits := 



P 



p-l 



log 



+ 2 



Moreover, from (jlSp we propose the following stopping criterion, instead of (jl3p : 



e-n-i 



< 10 



-r, (p-l)/p2 



(19) 



(20) 



5 Extrapolated Computational Local Order of Convergence (ECLOC) 

We give a relationship between A„ and p, a new technique to update the number of significant 
digits in an adaptive multi-precision arithmetic and a new stopping criterion. 

We start by deriving an expression of ECLOC as a function of the local order of convergence. 

Proposition 5.1 // A„ is the ECLOC defined in ^ and p > (j) is the order of convergence, then 

log\C\ 



Xr 



P{ 1 + 



(2p - 1) log \en-2\ 



21 



where C is given in 



Proof. Now, we write An in terms of en-2- To do that, we express log \en\ as a function of Cn-i 
and en-2- 



log I en I 



We obtain 



An 



2 log I en - e„_i I - log I e„ - 2en-i + e„_2 | 
2 log I en-i I + 0(1 e„/e„_i |) - log | e„_2 | + 0{\ en-i/e„_2 
log I e^_i/en-2 I + 0(1 en-i/en-2 I)- 

log I in I 



(22) 



log 


en-i/en-2 


+ 0(| en-l/en-2 1) 


log 


en-2/en-3 


+ 0(1 e„_2/en--3 1) 



From ([3]) we deduce e'^_i = (l + 0(e^_2)) , and taking into account (jlOp . 

C'^'''(in-2 (1 + Oie^Jj:^)). Next, we get 



log 



C^'4'-V(l + 0«-2)) 



log 



CVPetV''(l + 0«/-''2)) 



and the proof is complete. 
Proposition 5.2 Civen 



{2p - 1) log |en,2| + 21og |C| + 0(e^_2) 
(2p-l)log|en-2|+log|0| + 0(e^/_''2) 

— a and Cn = — oin, then 



□ 



^.n - C ^ e/ / ^'^P-^\ where /3 



2/9-1' 



(23) 
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Proof. Taking into account e„,_2 = C ^n-i + ^(^n-i)) ' write e„ in terms of e^-i: 



2er 



Ci/p,(2_P-i)/P(i + o(,./P). 



(24) 



where r = min{/9 — 1,0"}. Now, from (fM|) and e„_i = C ^^''ei'''' ^1 + 0(en'''')^ , we get 



(ji/p r^-i/pg Vp ('i + 0(e^/^) 



(2p-l)/p 



C(i-P)/P%M/P^ (^l+0(e„^/''')). 



(25) 



From (j25p . we have e 



{2p-l)/p2 



C '•^ ^^/^^ from which the proof immediately follows. □ 



Notice that ()23p allows us to implement an iterative method ([T]) with a multi-precision adaptive 
arithmetic. Instead of ([9]) we now consider the expression: 

n3 



Digits 



P 



X (-log| e„| + 2) 



2/9-1 

In addition, as an alternative to (|13p . ()23p provides the following stopping criterion 

leJ < l0-^(2p-l)/p^ 



(26) 



(27) 



6 Petkovic Computational Local Order of Convergence (PCLOC) 

In this section we provide a relationship between A„ and p. In addition, we derive a new technique 
to update the number of significant digits in an adaptive multi-precision arithmetic and a new 
stopping criterion. 

Proposition 6.1 If Xn is the PCLOC defined in ^ and p is the order of convergence, then 

'iog|cri-^l 



An 



P[ 1+ o 



P log \en- 



where C is given in ^ and T = f'{a). 
Proof. Setting 



f{xk) = r (efc + 0{el)) = Tck (1 + O(efc)) , k = n-l,n, 

and from ([3]), en = C'e^^-i (l + 0(e$^_x)i have 

log|/(x„)| plog|en_i|+log|rC|+0(|e-_i|) 



(28) 



(29) 



Xn 



log I /(xn-i)l log |e„_i| + log |r| + 0{\en-i\ 



Dividing the numerator and denominator in the right side of the preceding equation by log | e„_i| 
the proof is complete. □ 



Proposition 6.2 // we denote Qn 



fjXn) 
f{Xn-l)' 



then 



(30) 
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and (jlUp we have 



and the proof 

□ 



Proof. Taking into account 
immediately follows. 

The result (j30p allows us to substitute the error in Q by an expression that does not involve the 
exact root. Indeed, we implement the following adaptive multi-precision arithmetic scheme: 

^2 / /(x. 



Digits : = 



p-l 



log 



n— 1 ; 



+ 2 



Moreover, from ()30p we propose the following stopping criterion, instead of (jl3p : 



f{Xn) 



< 10 



'V (p-i)/p 



f{Xn-l) 

Table 1: Test functions, their roots and the initial points considered 





fix) 


a 


Xq 


{x^i , Xo} 




= - + a; - 2 


2.893289196304497788906356 


2.5 


{2.25,2.60} 




— x'^ + cos X — 2 


1.172577964753970012673333 


1.5 


{1.50,2.50} 




— 2 sin X + 1 ~ X 


2.380061273139339017212548 


2.5 


{1.00,2.00} 




= (x + l)e"-i - 1 


0.557145598997611416858672 


1.0 


{0.00,0.75} 






3.0 


2.94 


{2.90,3.10} 




— + cos a::. 


1.746139530408012417650703 


1.5 


{1.60, 1.90} 




= a:: — 3 In j: 


1 .857183860207835336456981 


2.0 


{1.00,2.00} 



(31) 



(32) 



7 Iterative methods and numerical results 

We consider in this section six iterative methods, c^^, A; = 1-^6, with local convergence order equal 
to 2, 3, 4, , 1 + \/2 and 1 + y/3 respectively. The first three methods are one point iterative 
methods without memory known as Newton's method, Chebyshev's method |11] and Schroder's 
method [12]. The other three schemes are iterative methods with memory, namely the Secant 
method and two variants (see [l3]). They are defined as 

(pl{Xn) = Xn- u{Xn), (33) 

<p2{Xn) = (j)l{Xn) - ^L{Xn)u{Xn), (34) 

4>3{Xn) = 4>2{Xn) - (^^L{Xnf - M{Xn)^ u{Xn), (35) 

(piiXn) = Xn - [Xn~l,Xn]j^ fiXn), (36) 

4'5iXn) = <p4.iXn) - [XnAiiXn]]]^ /{(l^^iXn)), (37) 

(peiXn) = (l)4{Xn) - [Xn,'^4>4{Xn) - Xn]]'^ fifpiiXn))- (38) 

where 

^(^) = ;^' ^(^) = 77^«(^X Mix) = ^^^uixf and i^'V]]' = f^y]~_"f^^y 

We point out that CLOC and PCLOC can be computed if n > 1, ACLOC if n > 2 and ECLOC 
if n > 3. If we have a method of higher order of convergence then multi-precision arithmetic 
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is required and will be used where necessary to obtain many correct figures. In general, from 
guess points X-j, . . . , x_i, Xq we obtain the admissible points Xi, . . . , Xj. Notice that if we use an 
arithmetic with at most 77 decimal digits, with stopping criterion |e/+i| < 10~^, then X7+1 will not 
be considered since it would need a mantissa with higher length to hold all the correct decimals. 
Hence, as x/ is the best admissible point, we will take n = / in the definitions of CLOC, ACLOC, 
ECLOC and PCLOC: 

Definition 5 The computational values considered in all numerical experiments are 

- log|e,| - logical ~ log|e^| - log|/(xj)| 

A = . 1 , A = — — 1 , A = pz — I and A = - — — . (39) 

log|ej_i| log|e^_i| log|ej_i| log|/(x^_i)| 

We have tested the preceding methods on seven functions using the Maple computer algebra 
system. We have computed the root of each function starting from the same initial approximation 
xo in (j33|) - ([35]) . and {x_i , xq} in ([36]) - ([38|) . Depending on the computational order of convergence 
used, CLOC (g]), ACLOC (P, ECLOC ^ or PCLOC jl]), we stop the iterative method when 
condition ([HI), ([20]), ([271) or ([32]) is fulfilled. Note that in all cases r/ = 2200. 

The set of test functions presented here were previously considered in [l3|. Table 1 shows the 
expression of each of these functions, their root with twenty five significant digits and their initial 
approximation. The latter is the same for all the iterative methods considered, considering one or 
two guess points depending on wether the algorithm works without or with memory. 

For each method and function, we have applied the four techniques with adaptive multi-precision 
arithmetic ([9]) , (fT9l) , (f26]l and ([3TI) . The number of necessary iterations to get the desired precision 
and the values of iterated points xi, . . . , x/ are the same. Table 2 shows the number of iterations 
needed to compute the root. In addition, the last four columns show the interval with minimum 
and maximum error produced in the computation of the corresponding Computational Local 
Orders of Convergence (CLOC, ECLOC, ACLOC or PCLOC) for the seven test functions. For 
instance, considering the CLOC and Newton's method (pi, let us denote X(A) the resulting error 
interval obtained in the computation of the CLOC: [min |Afc — p|, max |Afc — for each function 
/fc, fc = 1, . . . 7. 

Table 2: Min- max interval for error bounds 





/l 


h 


h 


h 


h 


h 


h 




J (A) 






T{\) 






T(\) 






TiX) 






12 


11 


10 


11 


12 


10 


11 


[2.8e- 


-5,l.le- 


-3] 


[3.7e- 


-5,1.5e- 


-3] 


[5.6e- 


-5,2.2e- 


-3] 


[6.0e- 


-5,1.2e- 


-3] 




8 


7 


6 


7 


8 


6 


7 


[8.9e- 


-5, 3.3e- 


-3] 


[l.Qe- 


-4,6.Qe- 


-3] 


[1.4e- 


-4,9.9e- 


-3] 


[2.1e- 


-4,4.5e- 


-3] 


<t>3 


6 


6 


5 


6 


6 


5 


5 


[8.8e- 


-6,1.3e- 


-2] 


[2.0e- 


-5,3.0e- 


-2] 


[3.5e- 


-5,5.1e- 


-2] 


[1.6e- 


-3, 1.2e- 


-2] 


4>i 


17 


18 


16 


16 


18 


14 


16 


[8.1e^ 


-6,5.8e- 


-4] 


[7.9e- 


-6,6.8e- 


-4] 


[1.2e- 


-5,9.4e- 


-4] 


[3.2e- 


-5, 5.5e- 


-3] 




9 


9 


9 


8 


10 


7 


8 


[5.5e- 


-5,3.0e- 


-3] 


[l.Qe- 


-4,4.4e- 


-3] 


[1.3e- 


-4,7.2e- 


-3] 


[1.9e- 


-3, 3.3e- 


-3] 


<^6 


8 


8 


7 


7 


8 


6 


7 


[3.6e- 


-5,3.7e- 


-3] 


[1.4e- 


-4, 1.6e- 


-2] 


[1.6e- 


-4, 1.7e- 


-2] 


[2.8e- 


~4,4.5e- 


-3] 



From these numerical tests, we can conclude that the CLOC produces the best approximations of 
the theoretical order of convergence of an iterative method. However, the knowledge of the root 
is required. Conversely, as we can see in the definitions of ACLOC ([5]), ECLOC ([7]) and PCLOC 
dS]), these parameters have the advantage that they do not involve the expression of the root a. 
Actually, in real problems we want to approximate the root which it is not know in advance. For 
practical purposes (see table [2] and Propositions 14. 1^ 15.11 and 16. ip we recommend ECLOC since 
it presents the best approximation of the local order. Nevertheless, PCLOC is a good practical 
parameter in many cases because it requires less operations. 
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A Numerical computations 



In this section we show the results of the six iterative methods, (l>k,k = 1-^6, cited above apphed to the test 
functions that appear in table [T] In each table we can find the necessary iteration number 7 to get the required 
precision; the error AA = |p — A| made by the CLOC A; the error AA = |p — A| made by the ACLOC A; the error 
AA = 1/9 - A| made by the ECLOC A and the error AA = |p - A| made by the PCLOC A. 

In tables |3]|8] we give the details of the error intervals observed in the computations of CLOC, ACLOC, ECLOC 
and PCLOC shown in table [1 



A.l Newton method method- 

The iterative Newton method defined by 

( given xo 

1 (?!)i(x„) = 



n > 0, 



(40) 



is used for the set of test functions (see table The errors AA, AA, AA and AA obtained are shown in tabled 



Table 3: The errors AA, AA, AA and AA for Newton's method 





I 


AA 




AA 




AA 




AA 




h{x) 


12 


1.803 e 


-4 


3.607 e 


-4 


2.404 e- 


-4 


1.086 e 


-3 




11 


2.790 e 


-5 


5.580 e 


-5 


3.720 e 


-5 


8.504 e 


-4 


Mx) 


10 


7.143 e 


-4 


1.430 e 


-3 


9.526 e 


-4 


1.220 e 


-3 


fi{x) 


11 


2.723 e 


-4 


5.448 e 


-4 


3.632 e 


-4 


6.446 e 


~4 


Mx) 


12 


1.109 e 


-3 


2.215 e 


-3 


1.478 e 


-3 


4.018 e 


-4 


Mx) 


10 


1.040 e 


-3 


2.082 e 


-3 


1.387 e 


-3 


1.121 e 


-3 


Mx) 


11 


1.512 e 


-4 


3.025 e 


~4 


2.016 e 


-4 


6.032 e 


-5 



A. 2 Chebyshev method -cl)2 method- 

The iterative scheme of Chebyshev method is 

given xo , 

, l/n ~ Xn — u[Xn) , {'^^) 

<})2(Xn) = yn - ]^L{Xn) u{Xn) , 71 >0. 

where u{x) = f{x)/f'{x) and L{x) = {f"{x)/f'{x)) u{x). The numerical results obtained for the set of test functions 
are listed in table ID 

Table 4: The errors AA, AA, AA and AA for Chebyshev's method 





I 


AA 




AA 




AA 




AA 




Mx) 


8 


2.077 e 


-4 


6.233 e 


-4 


3.739 e- 


'4 


3.048 e 


~3 


Mx) 


7 


7.185 e 


-4 


2.154 e- 


-3 


1.293 e- 


-3 


2.148 e 


-3 


Mx) 


6 


1.949 e 


~3 


5.858 e- 


-3 


3.511 e- 


-3 


4.527 e 


-3 


Mx) 


7 


8.917e 


-5 


1.437 e- 


-4 


1.033e- 


4 


2.109 e 


-4 


Mx) 


8 


3.318 e 


-3 


9.921e- 


-3 


5.965 e- 


3 


8.876 e 


-4 


Mx) 


6 


3.107 e 


-3 


9.350 e- 


-3 


5.600e- 


3 


3.612 e 


-3 


Mx) 


7 


2.017e 


-4 


6.051 


-4 


3.630 e- 


4 


4.643 e 


-4 



10 



A. 3 Schroder method -03 method- 

The iterative Schroder method is defined by 
' given xo , 

yn — Xn u(^Xn^ , 

' z„ = ^ L{xn) u{x„) , (42) 

(p3{x„) = Zn - {L{Xn))'^ — M{X„)J u{Xn), n>0, 

fix) f"(x) f"'{x) 2 

where ulx) = — L{x) = — ^—^ u{x) and M(a;) = — — --^ (tJ-ix)) ■ The errors obtained in the approximations 

fix) fix) Slf'ix) 

of p by expressions (|4]), ((5]), (O and (|8]) for the set of test functions (see table [l| are given in detail in table [S] 



Table 5: The errors AA, AA, AA and AA for Schroder's method 





I 


AA 




AA 




AA 




AA 




Mx) 


6 


8.809 e- 


-6 


3.524 e 


-5 


2.014 e- 


'5 


1.218 e 


-2 


Mx) 


6 


1.261 e- 


-3 


5.039 e- 


-3 


2.881 e- 


3 


2.091 e 


-3 


Mx) 


5 


2.142 e- 


-3 


8.585e- 


-3 


4.900e- 


3 


6.033 e 


-3 


Mx) 


6 


2.635 e- 


-4 


1.054 e- 


-3 


6.024 e- 


4 


1.567 e 


-3 


Mx) 


6 


1.299 e- 


-2 


5.129 


-2 


2.952e- 


2 


2.814 e 


-3 


Mx) 


5 


3.778 e- 


-3 


1.517 e- 


-2 


8.650 e- 


3 


4.521 e 


-3 


Mx) 


5 


6.252e- 


-5 


2.501e- 


-4 


1.429e- 


4 


3.027 e 


~3 



A. 4 The Secant method -cl)4 method- 

The fourth iterative function considered 04 is the Secant method defined by 

(given x-i , xo , 
Mxn) = ^"--FT~^\ TT^ — Tfixn), n>0. 
fiXn) - fix„~l) 

The errors gotten in the approximated computation of p by expressions ((31,(0,(171) and (ISJ for the set of test 
functions (see table H]) in the secant method are displayed in table 16] 



Table 6: The errors AA, AA, AA and AA for Secant method 





I 


AA 




AA 




AA 




AA 




Mx) 


17 


9.045 e 


-5 


1.466 e 


-4 


1.064 e- 


-4 


5.448 e 


-4 


Mx) 


18 


8.112e 


-6 


1.159 e 


-5 


7.925e- 


-6 


2.223 e 


-4 


Mx) 


16 


3.777 e 


-4 


6.100 e 


-4 


4.396 e- 


-4 


6.448 e 


-4 


Mx) 


16 


1.090 e 


-4 


1.788 e 


-4 


1.321 e- 


-4 


2.588 e 


-4 


Mx) 


18 


5.817e 


-4 


9.408 e 


-4 


6.811 e- 


-4 


2.107e 


-4 


Mx) 


14 


5.110e 


~4 


8.333 e 


-4 


6.098 


~4 


5.510 e 


-3 


Mx) 


15 


8.050 e 


-5 


1.295 e 


-4 


9.285 


-5 


3.187 e 


-5 
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A. 5 First variant of the Secant method -05 method- 

The iterative method <j>5 is written by the following equations 

given x-i, xq , 

, Vn = 04(a;„), j-^^^ 
Mxn) = yu- -^T^ ^ — r/(yn), n>0. 

I fiyn) - f{Xn) 

We apply this method to the set of test functions (table [1} and the errors AA, AA, AA and AA obtained are shown 
in table [7] 

Table 7: The errors AA, AA, AA and AA for (jis method 





I 


AA 




AA 




AA 




AA 




fi{x) 


9 


3.573 e^- 


-4 


8.632 e 


-4 


5.448 e- 


-4 


2.152 e 


-3 


h{x) 


9 


S.SlTe-- 


-5 


1.419 e 


-4 


1.026 


-4 


1.715 e 


-3 


Mx) 


9 


8.826 e- 


'4 


2.135 e 


-3 


1.349 e- 


-3 


1.507 e 


-3 


Uix) 


8 


5.266 e- 


'4 


1.271 e 


-4 


7.997e- 


-4 


1.247 e 


-3 


Mx) 


10 


1.536e- 


-3 


3.702 e 


~3 


2.337 e- 


-3 


5.563 e 


-4 


feix) 


7 


3.033 e- 


-3 


7.230 e 


-3 


4.375e- 


-3 


3.269 e 


-3 


frix) 


8 


4.718e- 


'4 


1.136 e 


-3 


7.120 e- 


-4 


1.876 e 


-4 



A. 6 Second variant of the Secant method -(f)e method- 
The iterative method (pe is defined by the following iterative scheme 



given 
yn 

(peixn) 



X-l , Xo , 

= 04(k„), 

= yn 



2(y„ - Xn) 



(45) 



/(2y„ - x„) - fix„) 



/(yn), n>0. 



The errors AA, AA, AA and AA of each sequence obtained by this method for each function of the set (see table [TJ 
are presented in table [8] 



Table 


8: 


The errors 


AA, AA, AA and AA for 


method 




J 


AA 




AA 


AA 


AA 


fiix) 


8 


4.641 e 


-4 


1.173 e-3 


6.377 e-4 


2.721 e-3 


f2ix) 


8 


3.614 e 


-5 


1.570 e-4 


1.448 e-4 


1.385 e-3 


f3ix) 


7 


2.626 e 


-3 


7.252 e-3 


4.481 e-3 


4.493 e-3 


fiix) 


7 


7.782 e 


-4 


1.477 e-3 


3.595 e-4 


1.705 e-3 


f^ix) 


8 


3.405 e 


-3 


8.740 e-3 


4.913 e-3 


1.122 e-3 


feix) 


6 


3.655 e 


-3 


1.672 e-2 


1.580 e-2 


4.021 e-3 


frix) 


7 


6.321 e 


-4 


1.870 e-3 


1.262 e-3 


2.826 e-4 
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